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Abstract 

Quasi-one-dimensional Peierls systems with quantum and thermal lattice 
fluctuations can be modeled by a Dirac-type equation with a Gaussian- 
correlated off-diagonal disorder. A powerful new method gives the exact 
disorder-averaged Green function used to compute the optical conductivity. 
The strong subgap tail of the conductivity has a universal scaling form. The 
frequency and temperature dependence of the calculated spectrum agrees with 
experiments on KCP(Br) and trans-polyacetylene. 
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The interesting new physics exhibited by quasi-one-dimensional electronic materials con- 
tinues to expand as advances in synthetic chemistry and semiconductor fabrication tech- 
nology make new materials available for study. Many quasi-one-dimensional materials 
undergo a structural instability known as the Peierls or charge- density- wave instability . 
Each chain has a periodic lattice distortion with twice the Fermi wavevector, 2kF, resulting 
in a gap 2A in the electronic spectrum at the Fermi surface, suggesting that such a gap 
should be clearly visible in the optical absorption spectrum. However, over the past twenty 
years, optical absorption measurements on a wide range of materials have consistently 
shown that there is a broad tail below the gap. Moreover, as the temperature increases the 
spectrum broadens considerably. There is is no accepted quantitative theory of the shape 
of the spectrum or its temperature dependence. While considerable effort has been made to 
understand excitation energies of solitons and polarons in these compounds, there is also 
no accepted theory for their absorption spectrum. 

The zero-point and thermal lattice motions have a significant effect on the electronic 
properties of the Peierls state and are, we believe, key to understanding the subgap absorp- 
tion A recent luminescence and Raman study of a metal halogen mixed- valence (MX 
chain) compound produced results consistent with this view [|IU|. If the relevant phonon 



frequency at g = 2kF, u!2kp, is much smaller than the optical frequency, as is the case in the 
optical absorption in most Peierls materials, the quantum and thermal lattice fluctuations 
can be modeled P,pA| by a static, Gaussian-random, backscattering potential, ^(x), with zero 



mean. Further, if the phonon dispersion near q = 2kp is ignored, ^(x) is fully characterized 
by the disorder-averaged correlator {C,{x)^*{y)) = 'y6{x — y). In general, the dimensionless 
disorder parameter rj = '-y/hvp^-, where fi? is the Fermi velocity, has contributions from 
extrinsic disorder such as due to impurities, r/e, as well as from the intrinsic disorder due to 
lattice fluctuations, r/j. The latter is related to the dimensionless electron-phonon coupling 
constant A, the phonon frequency at 2kF, and temperature 

viT) = + ,,(r) = + coth (||:) . (1) 



In this Letter, we develop a novel and powerful method for computing the electronic prop- 
erties of one-dimensional Peierls semiconductors with static disorder. We use this method 
to calculate the real part of the optical conductivity. We find (1) a strong subgap tail in 
the optical absorption; (2) the remarkable fact that the subgap conductivity, when prop- 
erly scaled, follows a universal scaling curve independent of the disorder parameter rj; (3) 
good agreement between the computed conductivity and the experimental conductivity for 
KCP(Br) IP] and trans-polyacetylene 0; (4) a universal scaling curve for the experimental 
conductivity for KCP(Br) at different temperatures with rj having the temperature depen- 
dence given by Eq. (|l]). Figure 1 contrasts the relatively good agreement with experiment 
achieved by our calculation against that for a Lorentzian convolution of the rigid-lattice 
conductivity [|1^. 



We consider the standard continuum model ||TT|jr3|| of noninteracting electrons in one 
dimension with a Peierls ground state, which, without the random term, has been studied 
extensively (see p. 795 in [R|) to understand the soliton and polaron excitations in conducting 



polymers. Our system is described by a Dirac-type equation |[T^ for the wave functions 
ipi{x) and 'ip2{x) for electrons moving with vp to the right and to the left, respectively, in 
the interval —L < x < L with a complex, Gaussian-random, backscattering potential ^(x): 
/ 
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The above equations correspond to the case where the Fermi wavelength of the electron is 
incommensurate with the lattice period |^|. Previous work has succeeded in computing the 



density of states and localization length for this model [jl6[. Here we report the first exact 
calculation of the Green function . 

We use (C^), the retarded (advanced) 2x2 matrix Green function to compute the 
real part of the frequency- dependent conductivity for huj ^ fc^T |n,18]: 

/-oo r -1 

j+^(E + huo, E\x -y)- j^^{E + huj, E\x - y)\ , (3) 
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where A is the cross-sectional of a chain and 



f^{E', E\x -y) = {hvFf{i:i \a^G^{y, x\E')a^G^{x, y\E)\ ). 



(4) 



(T3 is a Pauh matrix. Here we report the computation of the conductivity using the approx- 
imation {GG) ~ {G){G), which we expect to be vahd for hu ^ j/hvp, i.e., for optical but 
not dc conductivities. However, our method can be used to compute (GG) directly. 

Exact calculation of the Green function. Two linearly independent wave functions ijj 
and satisfy the boundary conditions: iIj{—L) = 1 and ipiL) = 1. These ijj and ijj will be 
statistically independent in the limit L 00 [0]. For x > y, 
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where {1^11^2 — ^2V'i) is a constant independent of x. We obtain G^ (G^) by solving for 
wave functions with Imii^ a small positive (negative) number. Then the denominator can be 
Taylor-expanded and 
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Next, we take the L ^ 00 limit and utilize the statistical independence of ip and ip to 
factor the average of the Green function into products of averages of G (Z)) and y. The 
quantities G„ = {Cn{x,y)) and Dn = {t)n{x,y)) are functions of {x — y) and yn = {yn{x)) 
is independent of x. Thus the averaged Eq. (^ can be written in terms of ?/„, G„(a; — y) 
and Dn{x — y). The following equations for yn and G„ can be derived from Eq. (^ using 
standard Fokker-Planck methods |T^, subject to the conditions yo = 1 and G„(0) = yn, 
which follow from Eq. ([7p for y^, and from Eq. (|^ and the equality {{ip2/4'i)"') = ((V'i/'^2)") 

g for Gn{oy. 
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2^ + ir]n ) y„ - yn+i - Vn-i = 0, (8) 
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A d{x — y) 

Dn satisfies the same equation as C„ with the initial condition -Dn(O) = Un+i- 

The key to our method is the fact that C„ and Dn decay very rapidly and exponentially 
as a function of n. The decay rate is determined by the dimensionless disorder parameter 
Tj = '-f/hvp^'- the larger t], the more rapid the decay. To solve Eqs. (|^) and @ numerically, 
we truncate the number of equations by setting ?/„ = for n > A^. For example, at rj = 0.1, 
the Green function is converged to 1 ppm for = 30. The method can be generalized to 
the calculation of (GG). The number of resulting equations scales as A^^. 

We calculate the conductivity using the exact Green function obtained by solving Eqs. 
and and the approximation (GG) ~ {G){G). In Fig. 2(a), we show the result for several 
different rj values. The t] = curve shows no absorption below the energy gap and a 
divergent behavior just above the gap. For finite rj values, the singularity is absent, and the 
spectrum broadens rapidly as rj increases. In the inset of Fig. 2(b), we plot the dependence 
on rj of three quantities characterizing the shape of the conductivity curve: F, which is the 
half-width for the low frequency side of the peak, the inverse of the conductivity at the 
peak, and 6u!, the shift of the peak frequency from 2 A. Quite remarkably, when we scale 
the conductivity by the peak value and the frequency by F, we find that the scaled curves 



for all rj values have a universal form independent of rj below the peak frequency pO| , pT |. 

Experimental data on the platinum, linear-chain compound K2Pt(CN)4Bro.3-3Il20 or 
KCP(Br) [^] show a similar scaling behavior. In Fig. 3, we plot the data for five tempera- 
tures. The scaling curve below the conductivity peak has approximately the same form as 
the theoretical curve in Fig. 2(b). We deduce the effective disorder parameter for each tem- 
perature by choosing rj such that the computed F(?7)/ci;pcak(^) [Hi is equal to F(T)/ci;pcak(^)- 
This criterion was used to choose the value rj = 0.43 in Fig. 1. In the inset of Fig. 3 we 
plot the 7] values thus determined versus temperature. The temperature dependence of rj 



is compared with the postulated form (|I]): an extrinsic temperature-independent rje plus 
an intrinsic temperature-dependent r/j. The fit value of extrinsic disorder is large. This is 
not surprising since x-ray and neutron scattering studies have shown that the bromine ions 
randomly occupy two different sites |2^. This disorder could cause significant scattering of 
the electrons along the platinum chains |^ . 

We have analyzed the optical conductivity data for trans-polyacetylene (CH)^; at room 
temperature [0. Since the phonon frequency at 2k p for this material is at a much larger 
energy scale than the room temperature, the optical conductivity does not depend sensi- 
tively on temperature. Consequently, unlike in KCP(Br), we cannot determine the separate 
contributions from extrinsic and intrinsic disorder. In Fig. 3, we show that the scaled (CH)^. 
sub g ap conductivity curve agrees well with those of KCP(Br). 

The peak of the conductivity is the product of ctq = e^vp/A/S. and the dimensionless 
ratio o"pcak(^)/o"o determined from fitting the experimental curves. Accordingly the peak 
conductivity yields ctq and hence vp since the gap 2A and cross-sectional area A are known. 
Table I shows that the fi^'s obtained this way are smaller than estimates based on band 
structure calculations. The values of rji{T = 0) for KCP(Br) and ri(T = 0) for (CH)^; 
obtained from the fitting are compared with estimates from the lattice zero-point motion. 

In conclusion, we have developed a new method for computing the electronic properties 
of quasi-one-dimensional Peierls systems including disorder and lattice fiuctuation effects 
and used it to obtain the absorptive conductivity at all temperatures. Our results compare 
well with data on KCP(Br) and trans-polyacetylene. The subgap conductivity curve has a 
universal scaling form. 

We thank Leo Degiorgi for helpful discussions and showing us his results before publi- 
cation. This research has been supported in part by the US DOE-Basic Energy Sciences, 
Division of Materials Sciences and the OSU Center for Materials Research. 
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FIGURES 

FIG. 1. Comparison of the real part of the frequency-dependent conductivity, (t(w), obtained 
from our calculations for a disorder parameter r] = 0.43 with experiment on KCP(Br) at T = 40 K 
1^] and the conductivity obtained from a Lorentzian convolution of the rigid-lattice conductivity 
for a dimensionless damping parameter ^ = 0.42 [12|. The frequency and the conductivity are 



scaled by the frequency at the peak, Wpeak) and by cipeak = '^{'-^peak), respectively. 

FIG. 2. Universal form of the computed subgap optical conductivity, (a) Conductivity spec- 
trum for different values of the disorder parameter rj (Eq. (||)). The spectrum broadens as rj 
increases. The conductivity scale is set by do = e'^vp/AA, where A is the cross-sectional area of a 
chain, (b) Scaling plot of cr{uj), where the frequency is scaled by T{rj), the half-width for the low 
frequency side of the peak. Note we claim only that the subgap {lo < ujpcak) conductivity scales. 
Inset: Dependence on the disorder parameter r] of Upeak; F, and 5io{r]) = 2A — Wpeak(f/) • 

FIG. 3. Scaling plot of experimental data for KCP(Br) |^ and trans-{Cli)x |0 at different 
temperatures. For each temperature, the corresponding disorder parameter value is obtained by 
choosing r] such that the theoretical T {-q) / u!pca.kiv) value is equal to T (T) / u!pca.k{T) ■ Inset: Com- 
parison of the temperature dependence of rj for KCP(Br) with Eq. (|^). We fit the data with a 
curve in which the first number (0.3) characterizes extrinsic disorder effects and the second number 
(?7i(0) = 0.07) depends sensitively on the value of uj2kp (twice the third number in K) which we 
took from experiment [^. The possibility of the scaling curve extending to temperatures above the 
three-dimensional ordering transition at Tp = 120 K fsl is based on the persistence of a 'pseudogap' 



above Tp |25|. 
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TABLES 

TABLE L Fermi velocity vp and disorder parameter ij determined from fits of experimental 
optical conductivity to our calculations. The peak height determines vp which is compared to 
the results of band structure calculations ■ The ij values obtained from the fits are compared 
to the intrinsic r]i{T = 0) = {6u/uo)'^{Aa/hvF), calculated using the band structure vp, and the 
zero-point motion Su and lattice distortion Uq in Table I of Ref . |^] . The electron-phonon coupling 
A is calculated from the estimated value of r]i and Eq. (1). 
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